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1. Introduction 

We analyze within the statistical hadronization model the rapidity densities at mid-rapidity 
in the central Au-Au collisions at a/s nn =200 GeV and ^^=130 GeV measured by STAR and 
PHENIX collaborations iJUggBSSBSBB 13 & The analysis is carried out in the 



grand canonical ensemble with the supplementary y$ parameter. This version of the statistical 



hadronization model is described in detail elsewhere [13, 14]. 



Statistical model analyses at mid-rapidity at RHIC are often performed by fitting the chemical 



freeze-out parameters to the ratios of measured rapidity densities []15J, |16|, 17|, [18 „ |19p. It is easy 
to see that in a boost invariant scenario the ratios of rapidity densities represent the ratios of the 
multiplicities since 

N A _ J(dN A /dy)dy _ (dN A /dy)Jdy _ dN A /dy 
N B J(dN B /dy)dy (dN B /dy)Jdy dN B /dy 

in which the second step results from the fact that the rapidity densities, being boost invariant, are 
independent on rapidity. We adopt a different approach here and do not include any ratios of rapid- 
ity densities in the analysis as this can lead to a large bias in the resulting best fit parameters [f2C|]. 
Instead, as in our statistical model analyses at AGS and SPS [13, |4]], we fit a scaling parameter V', 



common for all particle species, characterizing the size of the hadron emitting source in the limited 
rapidity window around y=0. 

We have shown before [ |i"3| ] that performing statistical analysis in a limited rapidity window 
at SPS energies can artificially enhance the relative strangeness content in the system. This is due 
to the fact that heavier particles have more narrow rapidity distributions and thus a kinematic cut 
around mid-rapidity may lead to a situation in which relatively larger fraction of massive particles 
is taken into account. However, if a sufficiently wide Bjorken type of boost invariant region, much 
wider than the typical width of the pion rapidity distribution coming from a single cluster at kinetic 
freeze-out is formed around mid-rapidity, then the rapidity densities of different particle species 
are approximately constant in this region. In this case it might be possible to determine the charac- 
teristics of an average fireball at mid-rapidity without introducing a bias arising from the kinematic 



cut [14]. 



2. Analysis and results 

The rapidity densities in 0-5% most central Au-Au collisions measured by STAR at ^^=200 
GeV and statistical hadronization model predictions for the same rapidity densities along our best 
fit parameters are shown in Table [T[ All other rapidity densities are corrected for the weak decay 
contamination except proton and antiproton rapidity densities which include the weak decay prod- 
ucts of A and A. Fit quality is generally good and rapidity densities are well described within the 
statistical hadronization model. Our findings agree with the similar analysis performed recently by 



the STAR collaboration [10] 



The STAR data in Au-Au collisions at ^^=130 GeV needs to be somewhat manipulated, 
since some of the rapidity densities are measured in different centrality windows. In order to 
keep consistency with the STAR 200A GeV data and PHENIX 130A GeV data, we have chosen to 
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particle 

K 1 

P§] 
0I§ 



A [10] 



A [10] 



S- [10] 



s+ ni 



v / r 3 e -o.7/r 

T [MeV] 
Ms [MeV] 

75 

% 2 /dof 
(E)/(N) [GeV] 



dN/dy 
322±32 
327±33 
51.3±7.7 
49.5±7.4 
34.7±6.2 
26.7±4.0 
7.70±0.30±0.85 
16.7±0.2±1.1 
12.7±0.2±0.9 
2.17±0.06±0.19 
1.83±0.05±0.20 
0.53±0.04±0.04 



SHM 
325.0 
327.4 
57.1 
53.5 
42.9 
30.9 
7.10 
16.0 
12.1 
1.87 
1.53 
0.63 



12.5±0.7 
161.0±3.9 

30.0±9.8 
1.02±0.05 

12.5/8 
0.99±0.04 



Table 1: Top panel: STAR rapidity densities in Au-Au collisions at ^^=200 GeV in the 0-5% most 
central collisions. Contribution from weak decays have been subtracted for all particle species except that 
proton and antiproton rapidity densities include weak decay products of A and A. Errors quoted are statistical 
+ systematic and are taken into account in quadrature. For pions, kaons and nucleons, statistical errors are 
negligible and thus are not shown. Bottom panel: The statistical hadronization model best fit parameters 
at chemical freeze-out and the mean energy per particle, averaged over all particle species, determined from 
the above experimental data. 



extrapolate all rapidity densities into the same centrality selection, namely to the 0-5% most central 
collisions. This is done by assuming that all final hadronic rapidity densities scale linearly 

dN dN h - 

~r = a + b ~i — t 2 - 1 ) 

ay dr\ 

with the negative hadron pseudorapidity density dN^- /drj in central and semi-central collisions. 
This is approximately true in the boost invariant scenario around mid-rapidity provided that soft 
processes dominate the particle production and all particle multiplicities are proportional to the 
negative hadron multiplicity. This seems to hold true relatively well for all the measured particle 
species. 

Using the STAR inelastic cross section versus dN^-fdri data from [g], we have calculated 
the negative hadron pseudorapidity densities dA^-/dT] in several centrality windows in the central 
pseudorapidity region (|tj| < 0.5) in order to estimate the rapidity densities for each of the particle 
species in the most central collision window. Particularly, we have calculated dA^/j /dr] = 297 
in the 0-5% most central collisions. This is in good agreement with STAR collaboration estimate 
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ANh- /drj 300±6% H |j] and we have chosen the same 6% as our error leading to the dN^- /drj = 
297±18 in the 0-5% most central collisions. 

Pions and lambdas are already measured in the 0-5% most central collisions and need no 
manipulation while the rapidity densities for other particle species need to be extrapolated. The 
negative hadron pseudorapidity densities as a function of centrality are published along kaon ra- 
pidity densities ^ in various centrality windows and thus we have fitted Eq. ( |2.1| ) to the data and 
used dNf,- /drj = 297±18 to estimate the K rapidity densities in the 0-5% most central collisions. 
Protons and antiprotons ^ are measured in the same centrality windows as kaons and thus we 
have used the same negative hadron pseudorapidity densities as for kaons in order to estimate the 
p and p rapidity densities in the 0-5% most central collisions. 

S~ and E + [Q] as well as (j) [Q] are measured in three different centrality windows and sim- 
ilarly to kaons and nucleons, we have estimated rapidity densities in the 0-5% most central col- 
lisions by assuming linear dependence on the negative hadron pseudorapidity density. Unlike in 
the case of kaons and nucleons, the negative hadron pseudorapidity densities estimated by STAR 
collaboration in these windows deviate somewhat from the ones we have calculated. In order to 
keep the same overall scale, we have used our calculated pseudorapidity spectrum to estimate the 
Z ± and rapidity densities in the 0-5% most central collisions. Also in |^], the 0-10% most 
central rapidity density for Cl + CI is given and without better knowledge, we have estimated the 
Cl + Cl rapidity density in the 0-5% most central collisions by scaling the yields with a factor 
(dNf,- /dri)5% / (dNf,- /dri)io%=l.06. No weak decay corrections have been applied to any of the 
rapidity densities except 71* densities do not include weak decay products of A, A and K®. 

Table || shows the extrapolated (experimental) and the statistical hadronization model predic- 
tions for the rapidity densities in the 0-5% most central collisions as well as the statistical model 
best fit parameters in Au-Au collisions at ^^=130 GeV. Again, the fit quality is good and most 
of the rapidity densities are described well within the statistical hadronization model. The chemical 
freeze-out temperature is about the same and baryon chemical potential is few MeV higher than in 
central Au-Au collisions at 200A GeV. Unexpectedly the strangeness under-saturation parameter 
Ys seems to be significantly over unity. To cross check our results, we can analyze the rapidity 
densities in central Au-Au collisions at the same center-of-mass energy ^^=130 GeV measured 
by the PHENIX collaboration. PHENIX has measured a subset of the particle species measured by 
STAR, namely A^, p, p, A and A [ 12, [Tij |. Our extrapolated rapidity densities agree with the 
corresponding PHENIX values (see Tables || and [3]). No weak decay corrections have been applied 
to the PHENIX data which is the reason why 71* rapidity densities deviate among PHENIX and 
STAR data at the same beam energy. In all cases, strong and weak decays corresponding to the 
experimental definitions are taken carefully into account in our analyses. The statistical model best 
fit parameters estimated from the PHENIX data agree with the ones coming from a fit to the larger 
data sample measured by STAR. This rules out the possibility that our extrapolations with STAR 
data would bias the fit towards unexpectedly large 75. 

We have repeated the analysis with 75 fixed to unity both with the STAR and PHENIX rapidity 
densities at y / ? NN =130 GeV. It seems that the minimum of the % 2 distribution is relatively flat in the 
case of the PHENIX data and the rapidity densities can be described well (x 2 jdof = 2.0/3) with 
the statistical hadronization model in absolute chemical equilibrium as well. The other resulting 
chemical freeze-out parameters are virtually the same as in Table [| Setting y$= 1 with the STARra- 
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particle 


dN/dy 


SHM 


71+ |g 


239.0d=3.0±2.0d=10.0 


233.8 


l_J 


239.0±3.0±2.0±10.0 


236.6 


A® 


17.2Q±0.4±1.72 


17.98 


A[§ 


12.3±0.3±1.23 


12.5 


K+ 


45.75±0.6±6.0±2.85 


46.92 


K 


43.16±0.6±5.4±2.57 


43.36 


P 


26.15±0.23±5.8±1.63 


32.17 


P 


18.85±0.16±4.12dbl.l5 


21.96 


M 


2.13±0.14±0.20±0.11 


1.80 




1.78±0.12±0.17±0.12 


1.42 


+ A 


0.586±0.11±0.056±0.035 


0.702 




6.09±0.37±0.69±0.34 


6.75 


y/^-O.V/T 


8.46±0.45 




T [MeV] 


160.5±4.3 




Mb [MeV] 


35.6d=12.8 




7s 


1.20±0.08 




X 2 /dof 


7.0/8 




(E)/(N) [GeV] 


1.00±0.04 





Table 2: Top panel: Estimated rapidity densities in Au-Au collisions at ^^=130 GeV in the 0-5% 
most central collisions measured by STAR collaboration. Rapidity densities of pions and lambdas are mea- 
surements while rapidity densities of other particle species are our estimates. Errors for our estimates are 
statistical + systematic + our extrapolation and are taken into account in quadrature. Statistical and system- 
atic errors quoted are the experimental errors in the most central experimentally accessible window, which 
are 0-6% for K , p and p, 0-10% for S and H + H and 0-11% for 0. No weak decay corrections have 
been applied except rapidity densities do not include weak decay products of A, A and K { }. Bottom 
panel: The statistical hadronization model best fit parameters at chemical freeze-out and the mean energy 
per particle determined from the above experimental data. 



pidity densities leads to somewhat larger temperature T = 169.1 d=4.4 MeV while baryon chemical 
potential is essentially unaffected. The anomalous behavior of strangeness phase-space occupancy 
factor at y / i NN =130 GeV is not yet well understood and the issue deserves further consideration. 

3. Discussion and Conclusions 

The resulting chemical freeze-out parameters can be compared with the ones we have deter- 



mined in central heavy-ion collisions at lower beam energies at AGS and SPS []14|]. One should 
notice that the previous analyses have been performed with full phase-space multiplicities while the 
statistical hadronization model analysis at RHIC employs rapidity densities measured in a limited 
rapidity window around mid-rapidity. Nevertheless, statistical hadronization model parameters in 
central Au-Au collisions at RHIC seem to be in good agreement with our interpolating curves for 
statistical hadronization model parameters determined at lower beam energies, see Figure 1. The 
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particle 


dN/dy 


SHM 


K+ 


276±3±36 




264.4 


K~ 


270±3.5±35 




269.6 


K+ 


46.7±1.5±7.0 




46.2 


K 


40.5±2.3±6.1 




42.9 


P 


28.7±0.9±4.0 




29.6 


P 


20.1±1.0±2.8 




20.6 


A 


17.3±1.8±2.8 




15.9 


A 


12.7±1.8±2.0 




11.8 


y> T 3 e -0J/T 


8.07±0.11 




T [MeV] 


158.0±5.9 




ll B [MeV] 


33.5±17.8 




Ys 


1.24±0.22 




X 2 /dof 


0.5/4 




(E)/(N) [GeV] 


0.97±0.06 





Table 3: Top panel: Rapidity densities around \y\ = in central Au-Au collisions at V / ^ NN =130 GeV 
measured by the PHENIX collaboration. All data taken from [O] except lambdas from [ 1 1 1. No weak decay 



corrections have been applied to any of the particle species. Errors quoted are statistical + systematic and are 
taken into account in quadrature in the analysis. Bottom panel: The statistical hadronization model best fit 
parameters at chemical freeze-out and the mean energy per particle determined from the above experimental 
data. 



solid line indicates the beam energy region in which the interpolating curves are determined while 
the dashed part is an extrapolation to RHIC energies. For the explicit formulas and details, see [p4|]. 

Figure 2 shows the strangeness phase-space occupancy factor y$ as a function of the center-of- 
mass energy of a colliding nucleon pair in central heavy-ion collisions. One can see that strangeness 
seems to reach absolute chemical equilibrium in central heavy-ion collisions at RHIC energies only. 
From the right panel of Figure 2 one can see that the RHIC systems at 200 and 130A GeV seem to 



follow the chemical freeze-out curve (E) / (N) = 1 GeV [21]. 



References 

[1] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 92 (2004) 1 12301. 

[2] J. Adams et al. [STAR Collaboration], Phys. Lett. B 612 (2005) 181. 

[3] C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 87 (2001) 1 12303. 

[4] C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 89 (2002) 092301. 

[5] C. Adler et al. [STAR Collaboration], Phys. Lett. B 595 (2004) 143. 

[6] C. Adler et al. [STAR Collaboration], Phys. Rev. C 65 (2002) 041901. 

[7] J. Adams et al. [STAR Collaboration], Phys. Rev. C 70 (2004) 041901. 

[8] J. Adams et al. [STAR Collaboration], arXiv:nucl-ex/031 1017. 



6 



Statistical hadronization model description for.. 



Jaakko Manninen 



180 

170 

160 

_ 150 - 
> 

g 140 r 
H 130 
120 
110 f- 
100 



1 



SPS Pb-Pb 




AGS 
Au-Au 



RHIC 
Au-Au 



10 100 
V s NN [GeV] 



700 
600 
500 
> 400 



m 300 



200 
100 








■ AGS Vr 


- 


Au-Au T 

\ 








\ 






RHIC 




Au-Au - 


SPS Pb-Pb 





10 100 
V Snn [GeV] 



Figure 1: Chemical freeze-out temperature (LEFT) and baryon chemical potential (RIGHT) as a function 
of the center-of-mass energy of a colliding nucleon pair in various central collision systems. The lines are 
our interpolating curves for central heavy-ion collisions [|l4||. 
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Figure 2: LEFT: Strangeness under-saturation parameter 7s as a function of the center-of-mass energy of a 
colliding nucleon pair in various central collision systems. 

RIGHT: The chemical freeze-out curve in central heavy-ion collisions. The line represents the condition 
(E)/(N) =1 GeV. 
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